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Abstract 

Superpositions of paraxial laser beam modes to generate atom-optical lenses based on the optical 
dipole force are investigated theoretically. Thin, wide, parabolic, cylindrical and circular atom 
lenses with numerical apertures much greater than those reported in the literature to date can be 
synthesized. This superposition approach promises to make high quality atom beam imaging and 
nano-deposition feasible. 
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I. INTRODUCTION 



The field of atom-optics offers considerable potential in applied and fundamental physics, 
both for atom beam lithography (to create nano-structures) and for atom beam mi- 



croscopy [3|, 14] . Here, the use of the optical dipole force using far detuned laser light for the 
manipulation of neutral atoms is considered. In this regime it yields a conservativepotential 
for the manipulation of atoms that is proportional to the laser light intensity l|, [sj. 

Already in 1978 Ashkin's group demonstrated neutral atom beam focusing using the 



optical dipole force [6]. Many techniques to focus atomic beams have been tried since: 
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the optical dipole force [1, 

Amongst optical dipole force approaches there are schemes based on pulsed laser config- 



urations [21 



221 ]. light confined by nano-apertures |18j, single-mode hollow beams |6l. l23l. l24 ] 
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or standing wave setups that yield tightly spaced ridges of the atomic deposition pat- 
terns 
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281 ]. Standing wave pattern approaches can also yield other deposition 



301 ]. but because of the high spatial frequencies involved, smooth profiles such 



32] 



as those desired for aberration-free atom-lenses wider than 200 nm cannot be synthesized 
with this approach [jj, M, 30, 31]. 

nn 

In the case of standing wave setups [1|, 130|] spherical aberrations give rise to pronounced 
pedestals, filling the gaps between patterned , |25|, |30| . This makes it impossibL 

lay down separate nano-wires. A pulsed approach should reduce the pedestal problem 
but remains constrained by the short spatial wavelengths typical for standing wave ap- 
proaches 33|]. A related approach 27|, that suffers less from pedestal problems, uses atomic 
de-excitation processes creating an effective transmission mask for excited noble-gas atoms 
to etch structures. Unfortunately, it appears to be unsuitable for direct deposition of metal 
atoms (they tend to stick to the deposition area regardless of their internal state). Its inher- 
ent filtering reduces atomic deposition rates and, more importantly, it does not redirect the 
center of mass of atomic motion and thus cannot be used for traditional imaging of atomic 
beams. 

Similar problems occur in the application of single-mode hollow laser beams 23|, l2J] as 
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optical imaging elements. Their waist is potentially wide, but their elongation leads to 
thick lenses with small numerical apertures: for realistic setups a diameter of the transverse 



parabolic part of the potential o: 



in the micrometer range 
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less than 200 nm arises in conjunction with focal lengths 



241 ] yielding unsatisfactorily small numerical apertures for 



atomic focussing. This implies that one would have to start out with already well focussed 
atomic beams and, yet, the resulting atomic point-spread function remains unsuitably wide. 
None of the approaches mentioned so far has been adopted as a solution for the problem 
of imagingof atomic beams in atomic microscopy [s], Ij] or direct atom-deposition litho- 
graphy ^i^ : a viable atom-optical lens still needs to be found. 




atom beam 



FIG. 1: (Color online) Arrangement of beams. 



Here, it is shown that only the superposition of many laser modes [34, l35|, |36| will allow 
us to generate wide atom-optical lenses based on the optical dipole force. We will find that 
widening the beams' waists is not a solution if atomic lenses with large numerical apertures 
are desired, because prohibitive increases in laser power are necessary. The idea of this paper 



is to superpose several odd Hermite-Gaussian TEMmn-modes, V^m,n [37|, l38|, l39| , such that 
all non-linear terms in the dependence of the electric field on the (transverse) x-direction are 
optimally suppressed, see Fig. [H This generates an electric field profile that varies linearly 
across a large part of the laser beam's cross section, see Fig. [21 and yields the desired 
parabolic laser intensity profile to generate an aberration-free atom-optical lens. 

After an introduction of the underlying idea in Section [Tll its possible implementation 
using Hermite-Gaussian modes to generate cylindrical lenses is described in Section UTTl Sec- 
tion |IV] generalizes this approach to a crossed beam configuration that yields thin spherical 
lenses. We conclude in Section |Vl 
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FIG. 2: (Color online) Electric field profile, i?2j+i(2;i 0, 0), at focal cross-section of Her mite- Gauss 
beams comprising superpositions of up to 23rd order odd modes (i.e. 2 J + 1 = 1, 3, . . . , 23; x-axis 
in units of beam waist wqx-, total cross-sectional beam power normalized to unity, Rayleigh lengths 
kept constant, eow£/2 set to unity). 

II. SUPERPOSITIONS OF ODD MODES 

We now consider cylindrical atom-lenses with a parabolic modulation in the a;-direction, 
see Fig. [U most of what follows can be translated into the scenario of circular lenses for 
which atomic beams co-propagate with the focussing laser beams 0, on their optical 
axis - instead of crossing through it. Such circular lenses would require the use of Laguerre- 



Gaussian instead of Hermite-Gaussian modes 23|] but they have the disadvantage of yielding 



either tiny lenses (in the case of strongly focussed laser beams) or thick lenses (for less 
focussed laser beams) We therefore do not investigate setups with laser beams co- 

propagating with the atomic beam here; instead, we will show in Section IIVI how to create 
a thin spherical lens using a combination of two orthogonally crossed multimode Hermite- 
Gaussian laser beams. 

A. Hermite-Gaussian modes 



Let us consider modes, ipm,n |37|, l38|, l39[, with transverse beam coordinates x and y 



propagating in the 2;-direction. The Rayleigh lengths z^i^ and zr^ associated with the two 
transverse coordinates, x and ?/, can be different from each other (focussed by different 
cylindrical lenses, say). In this case two different associated beam waist radii, wqx and 
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Woy, and Gouy-phases, (f)x{z) and 4>y{z) arise. In the paraxial approximation the normahzed 
modes have the form 
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Here, r = (x, y, z) is the position vector, ujl the frequency of the monochromatic laser, 
= ojl/c = 2tt/Xl its wavenumber, and iPmi^) = HrrjiC) exp(— ^^/2)/ 2°^m\^ ^ (m = 
0, 1, 2, . . .), with the Hermite polynomials H„ 
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38, 



39| . The wave front radii R{z) 



+zj^) I z, the beam radii w{z) = wqa/I + z'^/zj^ with wo = ^/XlZr/tt, and the longitudinal 

391], (j){z) = axctanlz/ zr), are all parameterized by the beams' 



Gouy-phase shifts 
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Rayleigh lengths zr; strictly speaking by z^^ and zr^, respectively. 
In a configuration, such as that displayed in Fig. [H one can generate 
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36| a wide cylin- 



drical atom-lens using a laser beam with an electric field composed of a suitable combination 
of odd modes 

J 

" (2) 



j=0 

Here, the beam is modulated in the x-direction whereas for the ?/-direction the purely Gaus- 
sian lowest order mode (po is employed. Note that this allows us to make the lens 'thin' 
in the ^/-direction. With increasing cutoff, J, the superposition pattern becomes increas- 



ingly dephased due to the action of Gouy's phase 
Section lYl 

Following reference 



35|, this will be further investigated in 



the modes in Equation ([T]) yield an electric field which is polarized 
in the y-direction with a small contribution in the z-direction due to the tilt of wave fronts 
off the beam axis (x, y, z are the unit- vectors and 3? stands for real-part) 



(3) 



For beams that are not too tightly focused we neglect the transverse derivatives. The 
associated time- averaged light intensity distribution then has the form 38] 



/2j+i(r-) =eo {E2j+i{r,tf) 



ujI 



2J+1 



(r) 



(4) 
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B. Normalization and Intensity Scaling 



With the normahzed modes of Eq. ([T]) and assuming that the sum of the coefficients 
Yl |c2j+iP in Eq. ([2]) is unity we use the normahzation 

/oo /»00 2 /"OO /»oo 

/ dxdy\i!2j+i{x,y,z)\'^ = 2"/ / dxdyl2j+i{x,y,z)\'^ = 1. (5) 

Assuming vahdity of the Raman-Nath approximation of neghgible transverse motion of 
the atoms ((x, z) = const.) [40], the atoms experience the integrated intensity distribution 
of the laser field given by 



hj+i{x,z)= I dyl2j+i{r)-^^^^ ^ 



2 w^{z) 



(6) 



We note that this integrated intensity / of beams of fixed total power reduces inversely 

— 1/2 

proportionally to their width wqx, that is, their field amplitudes scale with w^J . Fur- 
thermore the field gradients diminish with Wq^. This implies that the effective curvature 
of the integrated laser light intensity, jV^E^p, responsible for atomic focussing scales with 
Wq^. We face an unfavourable cubic scaling with the beam width if we attempt to expand a 
laser beam transversally in order to widen the effective lens without weakening its refractive 
power. Additionally, as we will show below, pure modes have small useful areas to generate 
lenses, the combination of these two factors makes a pure mode approach unfeasible. It 
forces us to employ the mode superpositions studied here. 
C. Optical Dipole Force 

We assume that the interaction between atoms and the laser light is well described by 
a two-level scheme (excited state 'e' and ground state 'gf') in rotating wave approximation 
with effective atomic line width F and resonance frequency uj = Ue — ujg. This leads to the 
expression I{r) F^/ (2/5) = r2(r)^ for the Rabi-frequency as a function of the ratio of the 
local laser intensity /(r) and the transition's saturation intensity Is = 7r/icF/(3A^) [s], 
With sufficiently weak laser intensity I and sufficiently large detuning 6^^ = ul ~ uj of the 
laser frequency ujl from the atomic transition frequency 00, the AC-Stark shift gives rise to 
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a conservative optical dipole potential which, to first order in I /Is-, has the form 



hV"^ I(r) 

We can determine the atomic de Broglie wave number k of atoms with mass M and initial 
kinetic energy Kq = {hn^Y / {2M) in terms of their kinetic energy K. Disregarding Doppler 
detuning, and assuming the validity of the Raman-Nath approximation {Kq ^ f4;), this 
allows us to calculate the associated phase shift 



1^4>{x,z) « / dv(K{r) - K„) = :i^EEl I dvi^l^ - 1) (8) 




h y " ' V Kq ' lQ^/K~oIs6^ 
The dependence of Eq. (Q on the inverse kinetic energy implies that best performance is 
achieved for monochromatic atomic beams; the approximations are in accordance with the 
Raman-Nath assumption 40|. Work by Drewsen et al. 42j showed that, for an atom-lens, 
chromatic dispersion can be reduced by tilting the laser beam with respect to the passing 
atomic beam, but the focal plane would have to be tilted as well. Such a tilt, however, 
elliptically stretches out the atomic beam's point-spread function. 

Aside from spherical aberrations, there are detrimental noise sources due to spontaneous 
emission of photons and light fluctuations. These tend to increase with increasing laser 
intensity but can be decreased by increased detuning [Sj] or through the use of more compli- 
cated optical level schemes {mI- Further discussion of their influences is beyond the scope 
of this paper. 



III. CYLINDRICAL LENSES 



According to eqns. ([8]) and ([9]) parabolic optical potentials give rise to parabolic atom- 
optical phase masks, as is required for 'perfect' atom lenses. In other words, we want 
the y-integrated electric field profile to depend linearly on the x-direction, see Fig. [2l In 
order to achieve this we integrate out the y-component, see Eq. ([S]), then Taylor-expand the 
field profile and finally choose the coefficients in Eq. ([2]) so as to cancel terms non-linear 
in X. Using the first 2 J + 1 odd field modes all non-linear terms up to (2J + 1)^^ order 
can be cancelled. The determination of the coefficients, 02^+1, involves the solution of a 
linear equation system and is easily performed. For instance, for the third superposition 
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FIG. 3: (Color online) Amplitude coefficients C2j+i of Hermite-Gauss superpositions ^'2j+i of up 
to 33^^ order modes (2J + 1 = 1, 3, . . . , 33). 



field, \[^5, comprising Hermite- Gaussian modes ipi^, ip^fi and ip^^, the relative strengths of 
the coefficients are C3 = cil8-\/6/71 and C5 = Ci2-\/30/71. For a normalized superposition 
the coefficient ci should be chosen accordingly. The family of the first twelve normalized 
superpositions {\E'2j+i,2J + 1 = 1,3, .. . ,23} is displayed in Fig. [21 the associated set of 
amplitude coefficients C2j+i is shown in Fig. [3J 

Mode-superpositions extend the "useful" linear part of the field profile yielding wider 
parabolic intensity profiles. Figure H] demonstrates that the useful parabolic part in the 
focal intensity profile expands with the number of superposition modes 2 J + 1 according 



to the V2J + 1-scaling, expected for a harmonic oscillator 35|]. Note, however, that the 
refractive power of the wider lenses is reduced (its atom-optical focal length is lengthened), 
because wider lenses have reduced transverse field gradients, see Fig. |2] and discussion fol- 
lowing Eq. ([6]). In order to compensate for this loss of refractive power, we can increase the 
transverse field gradient through either laser beam focussing in the x-direction, or through 
an increase in laser beam power. Focussing in the y-direction makes no difference since only 
the integrated intensity / matters. In the next subsection we show how much the power has 
to be raised to keep the atomic lenses' refractive powers equal. Subsequently, in Section [TVl 
we will investigate focussing in the x-direction; we will see that Gouy-dephasing constrains 
this focussing, the lenses must not be shrunken below a certain limit. 
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FIG. 4: (Color online) Integrated focal intensity profiles /(x,0) of Hermite-Gauss superposition 
beams comprising up to 23rd order modes, compare Fig. [2] (same units as in Fig. [21 vertical bars 
mark location of position 0.57 • \/2J + 1 ■ wqx, confirming harmonic oscillator-scaling S^])- 




FIG. 5: (Color online) Integrated focal intensity profiles I{x, 0) of Hermite-Gauss beams comprising 
up to 23rd order modes, compare Fig. [2] (same units as in Fig. total beam power adjusted such 
that all profiles have same curvature at origin as the dotted line parabola). 

A. Increased Beam Powers Compensates for Lenses Widening 



If we increase the total beam power P2J+1 for wider beam profiles according to the ratios of 
the modes' transverse derivatives, P2J+1 = Pi\dx'^i{x,0,0)/dx'^2j+i{x,0,0)\'^, the weakened 
gradient is power-compensated for by increased laser power. This way all optical potentials 
give rise to atom-lenses with equal refractive power, see Fig. O The necessary beam power 
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FIG. 6: (Color online) The increase in total beam power needed to achieve the power compensation 
described in the text and displayed in Fig. [S] as a function of mode number (solid red line) scales 
approximately hke |y • (2 J + 1)'^/^ (dotted black line). 

increase to achieve this compensation is sketched in Fig. [Hi The power savings due to our 
multimode approach are quantified in Subsection IIII CI 

B. Decreased Rayleigh-Lengths Compensate for Lenses Widening 

Alternatively to the beam-power increases just discussed, we can keep the total beam 
power for all beams equal and shrink the higher-order superposition-beams' Rayleigh lengths 
through increased beam focussing in the x-direction. This also allows us to compensate for 
the gradient reduction observed in Fig. [2l The corresponding laser intensity profiles are 
displayed in Fig. [7] and lend themselves to an efficiency analysis of the invested laser power. 
The vertical bars in this figure mark the points, d2j+i, where each intensity curve deviates 
from the enveloping parabola (dotted line) by 0.74 percent. They delineate the useful areas 
of the lenses. Beyond a deviation of 0.74%, spherical aberrations distort the atomic point- 
spread function of an imaged atomic beam too severely. The filled-in areas under the curves 
in Fig. [7] represent the laser power fraction contributing to the atom lens in each case. 
Higher-order superpositions clearly allow us to use the laser power much more efficiently. 
This is quantified in the next subsection. 
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FIG. 7: (Color online) Integrated focal intensity profiles I{x, 0) of Hermite-Gauss beams comprising 
up to 23rd order modes, compare Fig. [21 and their 0.74%-deviation marks, d2j+i, which lie at 
relative positions d2j+i/di = 1.00,3.24,... ,8.70 from the origin, compare Table [J (same units as 
in Fig. [21 in contrast to Fig. [5] total beam power normalized to unity, but Rayleigh lengths zrx 
readjusted such that all higher-order superpositions match up with curvature of the first mode 
case ^'i, see text). 



C. Lens Quality and Power Savings 



The 0.74%-criterion was extracted from the work by Gallatin and Gould 23(] who con- 
sidered, for example, the use of a 0.1 Watt laser detuned by roughly 40,000 linewidths. To 
achieve acceptable performance, the effectively useful beam area was found to be only some 



2di =140 nm wide (for a laser beam with a 2wox = 2/im waist diameter [23|). In other words, 
pure laser modes yield only a small useful window (in order to fulfill the 0.74%-deviation 
criterion only about 2di/ (2wox) =140 nm/2yum 7% of a cylindrical lens diameter or only 
the central 0.49 percent area of a circular lens are useful). Most of the laser power is wasted 
in the wings if no suitable superpositions of higher-order modes are employed. In our case of 
a cylindrical lens based on the Hermite- Gaussian mode ipi^, very similarly, approximately 
di/wox =6% of the width of the beam is useful, see Fig. [3 Additionally to the quantification 
of the useful area of the lenses, see Table [J this waste is meaningfully quantified through 
the determination of the fraction of power £^2j+i the laser beam contributes to the 'useful' 
part of the lens profile. We define it as the ratio of the laser energy contributing to the area 
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TABLE I: Lens Parameters 1^2 j+i and <?2j+i; compare FiglT] 



2J + 1 1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 

d2j+iM 1-00 3.24 4.75 5.74 6.45 7.00 7.42 7.78 8.06 8.32 8.53 8.70 8.87 9.01 9.15 9.27 9.36 

£2j+i [%] 0.048 1.6 5.1 9.1 13 16 20 23 25 28 30 32 33 35 37 38 39 

£2j+i/£i 1 34 107 190 269 344 411 472 526 576 620 662 699 735 766 795 825 



between the deviation points —d2j+i < x < d2j+i, in terms of the total laser power, namely 

^ _ IZo dy /-I;;;, dx hj+ijx, y, 0) 

Figure [7] and Table [T] summarize and quantify our findings. Specifically, Table [T] allows us 
to compare values for a single-mode atom lens for which Si = 0.048% with the superposition 
approach. For example, compared to mode the relative power savings in case of super- 
position \E'33 is 825, this translates into a power utilization of £^33 = 0.048% x 825 = 39%. In 
general the details of this behaviour depend on the chosen quality criterion but the underly- 
ing scaling is straightforward to derive. The useful fraction of the laser beam is proportional 
to a linear integral over the intensity and therefore grows with the third power of the position 
of the deviation mark £^2j+i/£^i = ('^2j+i/'^i)^- 



IV. SPHERICAL LENSES 



We now want to investigate the constraints that arise when an identical copy of the laser 
beam that travels along the 2;-axis, see Fig. [T], is additionally sent along the x-axis such that 
their crossed configuration leads to the simultaneous application of two cylindrical lenses 
giving rise to the application of a spherical lens to the atomic beam. Either the laser beams 
a.e .lightly displaced along the .-a.i. o. they are sufficiently detpad from each other such 
that despite their spatial overlap no harmful interference occurs [43]. 

Gouy's phase (j){z) = aicta.n{z/zji) ^ z/zji, introduces relative phases between the modes 



within each beam. Since the Gouy- 



Dhase varies strongest near the beam focus we have to 



consider its mode dispersive effects 35|]. If the beam is very strongly focussed (small value 



of zr^) the dephasing away from the focus z = is so rapid that non-linear aberrations 
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2J + 1 2J + 1 



FIG. 8: (Color online) The left panel illustrates the behavior of the relative deviation of the 
intensity distribution A/ from zero as it approaches the 0.74%-deviation marks (top and bottom 
grid). Here, A/23 is shown for the crossed configuration of two laser beams travelling along z and 
X-axis respectively. The value of the Rayleigh length z/j^; at which we find that the oscillatory 
behaviour of A/ along a constant radial perimeter just exhausts the upper and lower limits set by 
the deviation marks allows us to determine the associated value of Zmin- The latter is plotted as 
a function of maximum mode number, in the middle panel (the filled in red area is the forbidden 
area of too tightly focussed beams). The values of Zmin(2J + 1) in turn determine the positions of 
the turning points 0.57 • \/2 J + 1 • wqx (top green line), the positions of the deviation-points (i2j+i 
(middle blue line) and the minimal beam widths woxC^J + 1) (bottom red line), depicted in the 
right panel in units of the laser's wavelength A^. 

degrade the desired linear field profile over the width of the atom-lens. In other words, 
a lower limit for the Rayleigh lengths 2;min(2J + 1) as a function of the number of used 
modes '2J + 1' has to be determined in order to guarantee moderate dephasing. Whereas 
the absolute values for this lower limit are hard to derive from first principles, we can still 
work out the correct scaling with the mode number: 

The electric field is proportional to the superposition of the modes including the Gouy- 
phase factors; this can be approximated by -E2J+1 oc Yl'jt^\ ^^ji'j^^^'^ ~ Yl^/=i\ '^ji^ji^ + 
ijz/zji^). The expansion coefficients are positive and the wave functions are real at the 
focus z = 0. Since the first order term is purely imaginary the integrated intensity has to 
depend on z quadratically: hj+iiz) = -^2j+i(0) ■ [1 + 4— -^2j+i + ^i^'^)]- The term F2J+1 
has a complicated dependence on the number of modes, but, containing the square of sums 
of the form Yl'j^^\ j^i'^j^ roughly proportional to (2J + 1)^. When we consider the 
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relative deviation of the intensity profile near the focus from the focal intensity distribution, 

the widths 
see Fig. HI 



A/ = ^^Jr^, we find A/aj+i oc 4^ ■ (2 J + l)^. Additionally, we know that the widths 



of the superpositions scale roughly like those of the harmonic oscillator 
namely z oc ^/2J + 1. For constant relative intensity deviations A/2J+1 this implies const. = 
^^i"*"^ ■ (2J + 1)^ or zr^ oc (2J + 1)'^/^. A numerical investigation, see Fig. [HI confirms 
^min(2J + 1) = 0.8 ■ \l ■ (2 J + 1)^/^ as a good estimate for a lower bound on zr^. This 
relationship has been checked numerically and holds for 15 < 2J + 1 < 55. There is no 
reason to believe deviations might occur for values of 2 J + 1 > 55, but for small values of J 
the assumptions used in the derivation of the scaling law do not hold accurately, see Fig. HI 
Instead, the expression Zmm(2J+ 1) = 10.5 ■ ■ (2J+ 1)^/^ gives a much better estimate for 
2;min(2J + 1) in the range of 1 < 2 J+ 1 < 13. These lower limits for zr^ imply that the beam 
focus is several wavelengths wide and a posteriori confirms that the paraxial approximations 
hold for all cases discussed here, since the largest beam opening angle conforming with the 
lower limits presented here turns out to be roughly 7.5° for superposition \E'3. 

V. CONCLUSIONS 

For a possible experimental implementation of the ideas presented here it should be 
emphasized that throughout the use of a repulsive (blue-detuned) optical potential has been 
assumed since it allows us to build focussing lenses with a dark center reducing detrimental 
spontaneous emission noise. Equivalent logic applies to 'concave' atomic lenses which would 
best be implemented in red detuning, with dark centers as well. 

The Raman-Nath assumption becomes progressively worse the larger the numerical aper- 
ture of a lens. Trajectory simulations show a 'downhill' drift of atomic paths that can partly 
be compensated for by slightly weakening the rise of the potential through the suitable sub- 
traction of higher-order terms that lead to slight non-harmonic modifications of the lens, 
improving its performance. Clearly, if such fine-tuning is considered, the approximations 
underlying Eq. (jlj) and Eq. ([7]) might not be permissible. These considerations are beyond 
the scope of this paper. 

The techniques for the coherent superposition of laser modes have been experimentally 



demonstrated, see e.g. references 3J, l36|] and citations therein. We have found here that 



using the mode-superposition approach allows for very considerable laser power savings and 
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lenses can be made wider than is possible with pure modes. We come to the conclusion that 
for the design of atomic lenses, based on the optical dipole force, it is possible and necessary 
to coherently superpose suitable laser modes in order to create wide thin parabolic lenses 
with large numerical apertures. 

Obviously the approach presented here can be applied for the manipulation of stationary 



atomic clouds just as well as for atomic beams 4^, |45 |. 
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